Isospin effects and the density dependence of the nuclear symmetry energy 



ON 

o 
o 

(N 



S.R. Souza 1 ' 2 , M.B. Tsang 3 , B.V. Carlson 4 , R. Donangelo 1 ' 5 , W.G. Lynch 3 , and A.W. Steiner 3 
1 Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
Cidade Universitdria, 
CP 68528, 21941-972, Rio de Janeiro, Brazil 
2 Instituto de Fisica, Universidade Federal do Rio Grande do Sul 
Av. Bento Gongalves 9500, CP 15051, 91501-910, Porto Alegre, Brazil 
3 Joint Institute for Nuclear Astrophysics, National Superconducting Cyclotron Laboratory, 
and the Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA 

4 Departamento de Fisica, 
Instituto Tecnologico de Aerondutica - CTA, 12228-900 
Sao Jose dos Campos, Brazil and 
5 Instituto de Fisica, Facultad de Ingenieria, Universidad de la Republica, 
Julio Herrera y Reissig 565, 11.300 Montevideo, Uruguay 
(Dated: July 10, 2009) 

The density dependence of the nuclear symmetry energy is inspected using the Statistical Mul- 
tifragmentation Model with Skyrme effective interactions. The model consistently considers the 
expansion of the fragments' volumes at finite temperature at the freeze-out stage. By selecting 
parameterizations of the Skyrme force that lead to very different equations of state for the symme- 
try energy, we investigate the sensitivity of different observables to the properties of the effective 
forces. Our results suggest that, in spite of being sensitive to the thermal dilation of the fragments' 
volumes, it is difficult to distinguish among the Skyrme forces from the isoscaling analysis. On the 
other hand, the isotopic distribution of the emitted fragments turns out to be very sensitive to the 
force employed in the calculation. 
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Investigations on the density dependence of the sym- 
metry energy in the nuclear multifragmentation process 
have been stimulated by the discovery [JQ that the ratio 
i?2i between the experimental yields Y(A,Z) of a frag- 
ment of mass and atomic numbers A and Z, respectively, 
produced in similar reactions with different isospin com- 
positions, henceforth labeled '1' and '2' follows a scaling 

law @, I, S, @, 0, I, 1 EI EL El, El [LI EE El . This 

isoscaling law provides parameters a and (3 that are de- 
termined from the property 



R21 = 



_ Y 2 (A,Z) 



Cexp(aN + /3Z) 



(1) 



Yi{A,Z) 

where C is a normalization constant [l[ . 

Studies of the density dependence of the symmetry en- 
ergy have used Eq. (fT]) to: (i) probe the dependence of 
statistical models on the symmetry energy p, H, 0, EH 
and (ii) to probe the isospin c omp osition of the region 
emitting the fragments [3|, [lj], [l8|. In this paper, we 
are concerned with the former issue. The motivation for 
statistical model studies arises from the relationship be- 
tween the isoscaling parameter a and the nuclear sym- 
metry energy, whose leading term at low temperatures 
was shown to be 15. 101 



a « 4C sym [{Z 1 /A 1 f - (Z 2 /A 2 ) 2 ] /T , (2) 

where C sym denotes the symmetry energy coefficient Q, 
T is the temperature of the system, Z% and As, i = 1,2, 



correspond, respectively, to the atomic and mass num- 
bers of the decaying source. 

To be useful, reactions '1' and '2' should be chosen 
to produce systems at approximately the same temper- 
ature and density. The ratio in Eq. ([1]) involves yields 
of the same isotope; distorting effects associated with 
the deexcitation of the primordial hot fragments may be 
similar in the two reactions and lead to an approximate 
cancellation in the ratio. Theoretical calculations sup- 
port this assumption for primary distributions calculated 
from equilibrium statistical models IE 14 , EH ■ There- 
fore, measurements of the isoscaling parameters, through 
fits based on Eq. (fT]), may provide valuable information 
on the symmetry energy. 

This assumption has extensivel y b een exploited in 
many works 
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which employed 

the Statistical Multifragmentation Model (SMM) [V^. In 
order to reproduce the measured a parameter, C sym has 
been appreciably reduced, compared with the values usu- 
ally adopted in this model. The main conclusion of those 
works 0, H, [E E3, EI 1 12, IH is that there seems to be 
an important decrease of the symmetry energy at low 
densities. 

Although this result is reasonable on physical grounds, 
other studies [1, 0, 0| also provided a sound explanation 
to this apparent reduction of C sym . They suggest that 
surface effects associated with the symmetry energy, not 
considered in Refs. 0, H, 0, E3, EH, E3, EH' may also 
lead to significant reduction of the a parameter. This 
explanation seems to be more reasonable since the model 



used in all these studies @, @, @, 0, 1, S, E3, EH, El, 
EH is based on binding energy formulae that evaluate 
the isoscaler volume, surface and Coulomb terms at the 
ground state (saturation) density. Although the volume 
occupied by the total system is much larger than that 
of the ground state source Vq, the lower density values 
used in these statistical calculations are due to the space 
between the fragments. Therefore, in this scenario, there 
should be no sensitivity to the density dependence of the 
bulk symmetry energy. 

In this work we investigate this issue in a consistent ap- 
proach in which the properties of the fragments are cal- 
culated through the Thomas-Fcrmi approximation (TFA) 
at finite temperature. This version of the SMM, named 
SMM-TF, was presented in Rcf. 0. In this way, the 
changes to the fragments' energies and densities at the 
freeze-out stage are consistently obtained in the frame- 
work of the TFA. Thus, the total Helmholtz free energy 
F of a given partition mode is written in the same form 
as in the Imp roved Statistical Multifragmentation Model 
(ISMM) [Hf, also used in this work: 



F(T) 



C C 



oul 



7 2 



(1+ X )V3 



F t , 



.CO 



(3) 



A.Z 



-Ba,z + flAT) 



where 



Ccoui Z 2 
(1 + X) 1/3 A 1 / 3 



Ftrans = ~T(M - 1) log (Vf /A|) + T log ( 5o ^ /2 ) (4) 



tY,N a ,z log (g A ,zA 3/2 



A,Z 



N A .z 



log(AU,zQ 



In the above Eqs., CcquL denotes the Coulomb coefficient 
of the mass formula [21|, Aq and Zq are the mass and 
atomic numbers of the decaying source, Ba,z represents 
the binding energy of the fragment, Na,z stands for its 
multiplicity and M — z Na,z- The freeze-out vol- 
ume V x = Vb(l + x) is kept fixed (x = 2) for all the 
values of the source's excitation energy E* . The fac- 

( 3 /2 " 

tor M — 1, rather than M, as well as Tlog [goA ' 



in the translational contribution to the free energy, arise 
from the subtraction of the center of mass motion from 
the partition function of the total system. The spin de- 
generacy factor is denoted by gA.z, = /m n T 
corresponds to the thermal wave-length, and m n is the 
nucleon mass. As in the ISMM, the free volume reads 



V f (T) = V Q (1 + X)-Y, N A,zV A , z (T) , 



(5) 



A,Z 



where Va,z{T) is the volume occupied by each fragment. 
However, in the SMM-TF, V A ,z{T) is given by the TFA 
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FIG. 1: (Color online) Symmetry energy coefficient of cold 
nuclear matter for the SLy4 and Gs Skyrme forces as a func- 
tion of the density p. For details, see the text. 



[201 ] and therefore it differs from the ground state value 
used in the ISMM. The internal Helmholtz free energy of 
the fragment f A z m the SMM-TF model is also given by 
this microscopic approach. These are the only two differ- 
ences between the SMM-TF and the ISMM, as discussed 
in Ref. [2(j. All the other ingredients are the same. 

O wing to the plethora of Skyrme forces in the litera- 
ture 22, 1 23f | , we do not perform a detailed study using all 
the existing parameterizations. Instead, we selected two 



of them which give distinct Equations of State (EOS's) 
for the symmetry energy, i. e., the Gs [3] and the SLy4 
[2f| forces. Figure Q] shows the density dependence of 
the symmetry energy coefficient of cold nuclear matter, 
SE sym /A = C sym S 2 , where 5 = (p n - p p )/(p„ + p v ) and 
Pn {Pp) is the neutron (proton) density. One sees that, 
although both forces agree for densities close to the satu- 
ration value (po), the differences at lower densities can be 
appreciably large. The Gs and SLy4 forces provide exam- 
ples of a strongly density dependent (stiff) and a weakly 
density dependent (soft) symmetry energy, respectively. 
Both have bulk isoscaler incompressibility moduli in the 
range of 230 to 250 MeV [H]. Therefore, they are well 
suited to the present study. 



We confine our attention to the 40 Ar, 40 Ca 



5S 



Ni sys- 



tems, for which measurements of the isoscaling parameter 
a have been reported recently [l3[ . The decaying sources 
considered in the calculations below correspond to 80% 
of the compound systems, as 20% of the matter is re- 
moved in order to take the preequilibrion emission into 
account. To make a connection to the work of Ref. [l3j |. 
we adopt in the following their assumption that the Z/A 
of the source is the same as the original system. 

We begin our discussion by examining whether the 
isoscaling property should still hold if the thermal ex- 
pansion of the fragments' volumes is taken into account. 
The microcanonical treatment employed in the SMM-TF 
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FIG. 2: (Color online) Isoscaling predicted by the SMM-TF 
using the SLy4 for E* /A = 6.0 MeV. In the upper panel 
the yields of the hot primary fragments are used whereas the 
lower panel displays the results after the deexcitation of these 
fragments. The lines correspond to the best fit of the results 
using Eq. {T}. For details, see the text. 



does not allow one to derive analytical expressions to in- 
vestigate this issue. Furthermore, the dependence of the 
free volume on the species multiplicities, Eq. lead 
to highly non-linear terms in the Helmholtz free energy 
rendering the traditional grand-canonical formulas [|| in- 
valid. By minimizing F with respect to the multiplicities, 
as is done in Ref . @ , one may nevertheless obtain formal 
expressions which suggest that the isoscaling property 
should still be observed in this case. 

This is indeed found in our numerical microcanonical 
SMM-TF calculations, as is illustrated in Fig. [5J which 
shows R-2i for the SLy4 force and E*/A = 6.0 MeV. Sim- 
ilar results are obtained for the Gs force and for other 
excitation energies. The magnitude of the corrections 
due to the deexcitation of the hot primary fragments can 
be estimated by comparing the upper and lower panels 
of Fig. [5] Here, we simulated the decay with a simplified 
Monte Carlo Weisskopf model, which includes the emis- 
sion of nuclei up to oxygen. The parameterization of the 
cross-section for the inverse reaction was taken from Ref. 
[2^| . Following Ref. [28| , we calculate the density of states 
from the entropies associated with j* A z . This provides a 
consistent link between the primary stage and the deexci- 
tati qn p rocess. In agreement with previous calculations 

[SEE 

our results also suggest that a is not strongly 
sensitive to the deexcitation of the primary fragments. 

Figure [3] shows the comparison between the a values 
obtained in the different SMM models used in this work 
and the available experimental data The results 

corresponding to the primary fragments (top panel of 
this picture) reveal that a is fairly sensitive to the ther- 
mal dilatation of the fragments' volumes as the behavior 
of a in the SMM-TF calculations is different from that 
given by the ISMM. Switching from the SLy4 to the Gs 
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FIG. 3: (Color online) Isoscaling parameter a as a function 
of the excitation energy of the source, for different Skyrme 
forces and the ISMM, before and after the deexcitation of the 
primary fragments. The data are taken from Ref. For 
details, see the text. 

force leads to small differences in a, primarily within the 
small energy range 7.0 MeV < E*/A < 8.5 MeV. This 
corresponds to the region where the average fragments' 
density is approximately 0.65po- For the ISMM, investi- 
gations of the connection between a and C sym at finite 
temperatures show that Eq. J2J) can be fairly inaccurate 
at high temperatures, although the main conclusion that 
both quantities are strongly correlated remains valid [2^ | . 

It is interesting to examine whether comparisons of 
the SMM-TF calculations to experimental fragmentation 
data can clarify questions regarding the adequacy of the 
model and the effective forces selected to describe the 
multifragment emission. The ISMM results seem to fol- 
low the experimental trends more closely than the SMM- 
TF calculations. This picture does not survive the deex- 
citation of the primary fragments. Although the changes 
to a are small, the values predicted by the ISMM are 
systematically lowered, so that after secondary decay it 
lies below the data and the SMM-TF. The differences be- 
tween the calculations with the SLy4 and the Gs forces 
are very small after the decay of the primary fragments 
and our results suggest that it is difficult to distinguish 
between them from the isoscaling analysis. 

Distinct Skyrme forces lead to different values of 
Va,z(T), according to their properties. This directly af- 
fects the Helmholtz free energy of the system through 
changes in Vf and indirectly through changes in the inter- 
nal free energies of the fragments. Therefore, it suggests 
that the fragment composition should be very sensitive to 
the thermal expansion of Va,z- This is indeed observed in 
Fig. [U which displays the isotopic distribution of selected 
primary hot nuclei for the breakup of the Ar + Ni sys- 
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FIG. 4: (Color online) Isotopic distribution of selected pri- 
mary fragments for the Ar + Ni system at E* /A — 8.5 MeV. 
For details see the text. 



tern at E*/A = 8.5 McV. It has been shown in Ref. |20j 
that, at low energies, the isotopic distributions predicted 
by the SMM-TF are narrower than those given by the 
ISMM. However, they are very different at higher excita- 
tion energies, as shown in Fig. 31 where, for this system, 
the peak of the distribution shifts toward the proton rich 
isotopes. The effect is enhanced for higher excitation en- 
ergies and it is more pronounced in the case of the SLy4 
force than for the Gs force. We have checked that, al- 
though these trends remain true, the differences between 
the Gs and the SLy4 forces are reduced after the de- 
cay of the primary fragments whereas they remain large 
enough to clearly distinguish between the ISMM and the 
SMM-TF calculations. However, in order to draw pre- 
cise conclusions, a deexcitation treatment that explicitly 
takes the feeding between known discrete states of these 
frag ments into account, such as that presented in Ref. 
2M, should be developed and applied to the SMM-TF. 

In summary, we have tested, self-consistcntly, the sen- 
sitivity of equilibrium mutifragmentation theories to the 
density dependence of the symmetry energy. Our statis- 
tical calculations, which consistently incorporate Skyrme 
effective interactions, suggest that the corresponding pre- 
dicted differences in the isoscaling parameter a are not 
large enough to allow one to distinguish between sym- 
metry energies with very different density dependencies. 
This observation does not impact the use of a to probe 
the asymmetry of the emitting system, a use that has 
been instrumental for investigations of the density de- 
pendent symmetry energy 0, IT?! llil ] . On the other hand, 
the isotopic distribution of the fragments produced in the 
reactions retains some sensitivity to the effective force 
employed in the calculations. We believe that constrain- 
ing the key SMM-TF inputs, i.e. the source composition 
and the excitation energy, through experimental values 
of additional observables besides the isoscaling param- 
eters, will be critical to precision comparisons aimed at 



constraining the symmetry energy with isotopic distribu- 
tions. 
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